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ABSTRACT 

A 

The  effects  of  the  noise,  the  dissimilarity,  and  the  prefilter- 
ing on  the  detection  of  the  P-pP  delay  time  and  the  successful  decomposition 
of  the  mixed  P-  and  -pP  wave  by  employing  the  cepstrum  analysis  are  studied 
here.  These  effects  have  been  investigated  by  extracting  the  direct  P-phases 
and  the  sui  face  reflected  pP-phases  from  several  presumed  underground  ex- 
plosions 1 eastern  Kazakh,  and  remixing  them  at  various  time  delays  with 
the  addi  ion  of  various  levels  of  the  realistic  seismic  noise.  The  quality  of 
the  cepstrum  analyzed  results  as  a function  of  the  signal-to-noise  ratio,  the 
delay  time,  and  the  dissimilarity  between  the  P-phase  and  the  pP-phase  have 
been  obtained  for  a number  of  such  events.  In  general,  it  is  found  that  the 
cepstrum  analysis  can  detect  the  P-pP  delay  time  as  short  as  0.  4 second  and 
successfully  recover  the  P-  and  the  pP-phases  for  the  signal-to-noise  ratio 
down  to  about  12  dB.  When  the  similarity  coefficient  between  the  P-  and  the 
pP-phase  5 is  less  than  1.0,  the  detection  of  the  P-pP  delay  time  becomes 
not  so  obvious  as  it  is  when  they  are  identical.  The  dissimilarity  between  the 
P-phase  and  the  pP-phase  will  affect  the  detection  of  the  P-pP  delay  time 
more  seriously  than  the  noise  in  the  cepstrum  analysis.  ^ Prefiltering  of  the 
noisy  mixed  P-  and  -pP  wave  with  the  standard  oandpass  filter  will  not  affect 
the  cepstrum  analyzed  results  as  long  as  the  filter  band  covers  the  entire  sig- 
nal spectrum;  although  it  seems  that  no  prefiltering  will  offer  better  detection 
of  the  P-pP  delay  time. 
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INTRODUCTION 

In  the  previous  reports  on  cepstrum  analysis  (Lane  and  Sun, 
1974a,  1974b,  1975),  the  noise  was  not  taken  into  consideration  explicitly, 
since  we  usually  chose  the  signals  with  reasonably  good  signal-to-noise  ratio 
in  order  not  to  complicate  the  problems  of  interest.  Surprisingly,  among  the 
few  published  articles  on  the  applications  of  the  cepstrum  analysis  to  seis- 
mology, there  was  only  one  occasion  where  the  noise  was  discussed  (Ulrych, 
1971).  Even  so,  it  was  only  given  as  an  example  using  randomly  generated 
white  noise  instead  of  the  real  seismic  noise  associated  with  the  signal,  and 
the  signal-to-noise  ratio  was  quite  fair  (20  dB). 

Just  as  in  most  signal  processing  problems,  we  can  expect  that 
noise  will  affect  the  capability  of  the  cepstrum  analysis  on  the  detection  of  the 
P-pP  delay  time  and  the  successful  decomposition  of  the  mixed  P-  and  -pP 
wave.  The  questions  to  be  asked  are  how,  and  to  what  degree,  is  the  effect 
of  the  noise.  The  answers  to  these  questions,  together  with  the  noise  level  at 
an  observational  array,  can  lead  to  more  meaningful  cepstrum  analysis  of  low 
signal-to-noise  ratio  signals  (which  are  usually  associated  with  the  low  magni- 
tude presumed  nuclear  explosions),  and  an  estimate  of  the  magnitude  and  depth 
at  which  these  events  might  be  concealed  by  overburying. 

In  addition  to  the  effec.  of  noise  on  the  cepstrum  analysis,  the 
effect  of  dissimilarity  between  the  P-phase  and  the  pP-phase  has  also  been  in- 
vestigated. The  information  obtained  from  the  comparison  between  the  two 
kinds  of  effects  can  be  used  to  explain  the  cepstrum  analyzed  results  for  the 
naturally  mixed  P waves  whci  1 the  pP-phase  usually  can  not  be  considered  as 
the  exact  replica  of  the  P-phase, 


1-1 


The  signal  and  the  noise  estimates  used  in  this  study  are  taken 
from  the  teleseismic  short-period  P waves  of  presumed  underground  explo- 
sions in  eastern  Ka/.akh  as  recorded  at  NORSAR.  Therefore,  it  will  not  be 
appropriate  to  generalize  the  experimental  results  in  this  report  for  some 
other  class  of  signals  and  noise  without  due  consideration. 

In  Section  II,  a theoretical  discussion  of  these  problems  will  be 
given.  The  experimental  results  using  the  synthetic  data  are  presented  in 
Section  III.  Finally,  the  results  from  the  analysis  of  several  naturally  mixed 
P-  and  -pP  waves  with  various  signal-to -noise  ratios  are  given  in  Section  IV. 
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SECTION  II 

THEORETICAL  CONSIDERATIONS 


In  this  section,  we  look  into  the  problem  from  the  mathematical 
point  of  view.  Although  the  derived  equations  for  the  complex  cepstrum  and 
the  resolved  signals  can  not  be  explicitly  expressed  as  a function  of  the  signal- 
to-noise  ratio  and  the  dissimilarity  of  the  two  signals,  we  can  see  from  the 
equations  where  the  noise  and  the  dissimilarity  enter  into  the  problem. 


A.  THE  EFFECT  OF  NOISE 

Let  us  consider  a mixed  signal  consisting  of  two  identical  sig- 
nals with  some  additive  noise  as  follows: 


x(n)  = s(n)  + as(n-nQ)  + «(n) 


(II-l) 


where  a is  a scale  constant,  nQ  is  the  time  delay,  and  »<n)  is  the  additive 
noise.  The  equation  (II- 1)  can  be  rearranged  in  the  following  form: 


x(n)  = s(n)  * m(n)  + f (n) 


(II- 2) 


, . c/  . . Sf„  „ \ ixfnl  is  the  Dirac  delta  function)  is  the  multi - 
where  m(n)  = 8(n)  + a8(n-nQ)  (6(n)  is  tne 

path  operator  which  convolves  with  the  signal. 

The  complex  cepstrum  of  x(n)  is  obtai  .ed  by  the  forward  cep- 

strum  transformation  as  given  below.  First,  take  the  z-tra.  sform  of  x(n) 


as  follows, 


X(z)  = S(z)  (1  + az  ) + E(z)  , 


(II- 3) 


where  S(z)  and  E(z)  are  the  z-transforms  of  the  signal  and  the  noise,  res- 
pectively. To  obtain  the  complex  cepstrum  of  x(n),  we  have  to  take  complex 
logarithm  of  equation  {II- 3)  as  follows: 


A 

X(z) 


”no 

£n  Js(z)l  + £n(l  + az  ) 


in 


1 + 


E(z)/S(z) 

-nQ 

1 + az 


(II-4) 


Now,  if 


max 

-n<(o<  IT 


-in  ai 

ae  ° <1  , 


the  complex  cepstrum  of  x(n)  can  be  expressed  in  the  following  form: 

k 


x(n)  = s(n)  + ^ ^ ( - 1) 


k+  1 a 


, 8(n  - kn  ) 

k o 


k-  1 

+ z-transform 


1 ffn(l  + F<*>'S[Z>  ) 

1_  ' 1 + az  n°'. 


(H-5) 


where: 


and 


z-transform’1^  J stands  for  inverse  z-transform 
s(n)  is  the  complex  cepstrum  of  the  signal. 


In  general  the  complex  cepstrum  of  the  noise,  which  is  represented  by  the 
last  term  in  equation  (II- 5)  can  be  expected  to  spread  over  the  entire  cepstrum 
time.  However,  with  good  signal-to-noise  ratio,  the  last  term  in  equation 
(II. 5)  can  be  expected  to  be  small  and  the  cepstral  peaks  (at  least  the  first 
few)  of  the  multipath  operator  (which  is  the  second  term  in  equation  (II-5)), 
should  be  readily  identifiable.  The  appropriate  comb  filter,  which  eliminates 
the  second  term  in  equation  (II-5),  can  be  applied  at  the  right  cepstrum  time, 
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i,  e.  , n = kn  . Then  the  resolved  signal,  which  is  the  inverse  cepstrum 
o 

transform  of  the  filtered  x(n),  will  be 

£(n)  = s(n)  + z-transform 

Therefore,  the  additive  noise  in  the  input  mixed  signal  x(n)  will  show  up  in 
the  resolved  signal  ^(n)  us  additive  noise,  which  is  the  original  raise  filter- 
ed by  the  linear  system  corresponding  to  (1  + az  n°) 

F rom  the  expressions  given  in  equations  ( II - 5)  and  (II-6),  it  is 
hard  to  see  explicitly  the  effect  of  the  noise  on  the  cepstrum  analysis.  That 
is,  the  explicit  dependence  of  the  detection  of  the  cepstral  peaks  and  the  quality 
of  the  resolved  signal  on  the  signal -to -noise  ratio  can  not  be  obtained  analy- 
tically. Nevertheless,  it  is  quite  clear  that  the  noise  will  definitely  smear 
the  cepstral  peaks  of  the  multipath  operator,  and  the  degree  of  the  smear 
should  depend  on  the  signal-to-r.oise  ratio. 


E(z) 


1 + az 


-nr 


(II -f>) 


B.  THE  EFFECT  OF  DISSIMILARITY 

The  dissimilarity  here  refers  to  the  situation  when  the  com- 
ponent signals  which  form  the  mixed  signal  are  not  identical.  We  may  treat 
this  case  in  the  following  noise-like  analysis.  Let  the  mixed  signal  be 


x(n)  = s^n)  + s^n  - nj 


(II-7) 


where  s (n)  and  s (n)  are  not  identical. 
1 £ 

analysis,  let 


Proceeding  as  if  this  were  a noise 


*d(n)  = s2(n)  _ asj(n) 


(II-8) 


where  a is  a constant 

II  2 
*d(n) 

1 

as  follows: 


and  can  be  uniquely  determined  by  requiring 
be  a minimum.  Then  equation  (II - 7)  can  be 


that  the 
rewritten 


II- 3 


x(n)  Sj(n)  + as ^(n  -nj  -t  tcj(n~n0)  • 


(IT-9) 


Comparing  equations  (II -9)  and  (II - 1 ) , it  is  plausible  to  treat  the  difference 
between  the  two  component  signals  as  the  noise.  We  must  keep  in  mind, 
though,  that  t (n)  is  closely  related  to  the  signals  s (n)  and  s (n).  In  gen- 

Cl  L L, 

oral,  this  is  not  the  case  for  the  noise  as  discussed  in  Subsection  A. 

Following  the  same  procedure  used  in  Subsection  A,  the  com- 
plex ccpstrum  of  x(n)  can  be  found  to  be: 


x(n)  = s 1 (n)  + (-l)k+1  y-  6(n-knQ) 


+ z-transform 


"M"  7 


Ed(z)z  °/51(z)' 


+ az 


(11-101 


where  E ,(z)  and  S,(z)  are  the  z-transforms  of  * (z)  and  s (n),  respec- 
d 1 d 1 

tively,  and  s^(n)  is  the  complex  ccpstrum  of  s^(n).  Equation  (II - 10)  is  sim- 
ilar to  equation  (II- 5),  but  now  the  last  term  in  the  equation  is  attributed  to 
the  fact  that  s (n)  and  s (n)  are  non -identical.  Again,  in  general,  we  can 
xpeot  that  the  last  term  will  be  non-zero  for  the  entire  ccpstrum  time  do- 
main. However,  when  s^(n)  and  s^n)  are  similar  |*d(n)|  can  be  con- 
sidered small  as  compared  to  js^(n)j  and  | s ^(n) | . Then  the  last  term  in 

equation  (II - 1 0)  can  be  expected  to  be  small  and  the  cepstral  peaks  of  the  multi- 
path  operator  might  be  identified.  The  application  of  the  appropriate  comb  fil- 
ter to  x(n)  w;  i yield  the  resolved  signal 

r e (z>  z'n°i 

#1  t \ , ^ . r - 1 I a / tt  in 


s.(n)  = s , (n)  + z-transform 

~1  1 


1 + az  n° 


The  estimate  of  the  s^n)  will  be 


(II-H) 


I 

B 


s(n)  = x(n)  - s (n) 
^ / 1 


.1 

- s (n  - n ) - ^-transform  _n 

2V  o 1 { az  0 


(11-12) 


Thus,  the  dissimilarity  of  S,H  and  s^n)  wil.  appear  in  the  cepstrum  re- 
solved  signals,  s](n)  and  s^n)  , as  .he  tillered  additive  noise. 

From  .he  above  discussion,  it  seems  that  the  noise  and  the  dis- 
similarity will  affect  the  cepstrum  analysis  in  the  s.m.lar  fashion.  In  practice 
however,  « ,<n)  in  equation  (11-8)  can  be  expected  .0  be  correlated  to  the  sig- 
nal more  than  the  noise  t(n)  and  the  spectral  content  of  *d<n)  will  e s.m 

1Ur  ,0  that  of  the  signal.  Hence,  the  complex  cepstrum  attributed  to  t <n>  , 
whi-t,  is  the  las.  term  in  equation  (11-10).  might  be  distributed  more  ,n  t e 
same  range  where  the  signal  lies  than  elsewhere  along  the  cepstrum  time  axis. 
Therefore.  is  conjectured  that  the  dissimilarity  might  affect  the  detection 
of  the  cepstral  peaks  of  the  multipath  operator  more  seriously  than  does  t e 


noise. 


C.  THE  EFFECT  OF  PREFILTERING 


For  the  noisy  mixed  signal  as  given  in  equation  (II-D.  in.uitive- 
ly  one  would  want  .0  filter  out  part  of  the  noise  in  the  signal  x(n)  .0  improve 
the  signal-to-noise  ratio.  But  the  tradeoff  in  improving  the  stgnal-to-no.se 
rati,  usually  is  the  signal  distortion.  In  the  cepstrum  analysis  the  signa 
distortion  might  make  it  harder  ,0  detect  the  cepstral  peaks  o,  the  mu, t, path 

ope  rato  r . 

The  spectrum  of  the  mixed  signal  x(n)  can  be  easily  obtained 

by  letting  z = ej“  in  the  equation  (II- 3) 


X(u>)  = S(w)  M(w)  + E(<o) 


II- 5 


where  10  is  the  angular  frequency  and  j = . S (to)  is  the  spectrum  of 

s(n)  and  M(cu)  is  the  spectrum  of  multipath  operator  m(n),  i.  e.  , 


-iwn  , 

M (co)  = (1  + ae  °) 


Writing  the  spectrum  in  the  polar  form 


i [df(uj)  + 1 

X(o»)  = | A(co)  | | M(oj)  | e *- 


(11-14) 


v'he  re 


A(o»)  = S(ai)  + . tf'M  is  the  Phase  of  A(w)* 

1 + ae"JtlJn° 


| M(  to ) | 


1 + a + 2a  coston0  , and 


-a  sinom 


1 + a coswn 


Thus,  in  both  the  amplitude  and  the  phase  spectra  of  x(n),  the  multipath  op- 
erator m(n)  represents  a periodic  modulation  whose  period  in  the  frequency 
domain  is  equal  to  the  inverse  of  the  delay  time  in  m(n).  Figure  11-1  shows 
the  amplitude  spectrum  of  x(n)  which  will  add  in  the  understanding  of  the 
periodic  modulation.  By  examining  equations  (11-13),  (11-4),  and  (11-5),  it  is 
not  hard  to  see  that  the  cepstral  peaks  of  the  multipath  operator,  which  is 
the  second  term  in  equation  (11-5),  are  the  direct  result  of  the  periodic  mod- 
ulation. Hence,  it  is  understood  that  those  ripples  of  the  periodic  modulation 
shown  in  Figure  11-1  carry  the  essential  information  for  the  detection  of  the 

cepstral  p^aks. 

Now,  let  us  apply  a standard  bandpass  filter  to  X(o»).  The 
filter  response  is 
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FIGURE  11-  \ 

amplitude  spectra  of  x(w),  a (u),  and  mm 


FIGURE  11-2 

FILTERED  AMPLITUDE  SPECTRA  OF  X(w) 
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AnM  = 1 


, 

«j  £ W < Oi 
L ri 

Ab(u,)  = 0 " < w > w„  ' 

The  filtered  X(o.)  is  shown  in  Figure  11-2.  It  is  clear  that  we  removed  the 
spectra  for  w < w and  <a  > , but  we  also  removed  the  ripples  in  those 

frequency  ranges.  Thus,  the  periodicity  of  the  amplitude  modulation  original- 
ly associated  with  M(ui)  is  destroyed.  This,  in  turn,  will  affect  the  nice  ap- 
pearance of  the  cepstral  peaks  as  predicted  by  the  second  term  in  equation 
(11-S).  Therefore,  it  is  reasonable  to  believe  that  any  prefiltering  technique 
which  will  destroy  the  periodicity  of  the  amplitude  modulation  produced  by 
the  multipath  operator  will  not  help  the  cepstrum  analysis  for  detection  of 

the  delay  time. 


SECTION  III 

EXPERIMENTAL  RESULTS  - SYNTHETIC  SIGNALS 


As  pointed  out  in  the  previous  section,  from  the  mathematical 
expressions  it  is  hard  to  see  explicitly  the  effects  of  the  noise,  the  dissimi- 
larity, and  the  prefiltering.  Therefore,  to  draw  some  practically  useful  con- 
clusions about  these  effects,  the  most  direct  way  is  to  cepstrum  analyze  a suf- 
ficient number  of  synthetic  signals  made  from  the  class  of  signals  in  which  we 
are  interested.  This  section  presents  the  experimental  results  to  explicitly 
demonstrate  these  effects  on  the  detection  of  the  P-pP  delay  time  and  the  suc- 
cessful decompositions  of  the  mixed  P-  and  -pP  wa'-es.  These  effects  are 
treated  separately  in  order  to  isolate  one  effect  from  the  other  for  better  de- 
monstration. The  class  of  signals  which  are  used  here  are  the  teleseisrmc 
short-period  P waves  (with  sampling  rate  of  10  samples  per  second)  of  the 
presumed  underground  explosions  in  the  eastern  Kazakh  (EKAZ)  recorded  at 
NORSAR. 

A.  PREPARATION  OF  SYNTHETIC  SIGNALS 

First,  four  high  signal-to-noise  ratio  EKAZ  events,  as  listed 
in  Table  III-  1 , are  cepstrum  analyzed  and  successfully  decomposed  into  their 
P-phases  and  pP-phases.  Figure  III- 1 shows  the  typical  result  for  KA7/081/ 
04N.  The  solid  line  is  the  resolved  P-phase  and  the  dashed  line  is  the  pP- 
phase.  These  resolved  P-phases  and  pP-phases,  together  with  the  real  seis- 
mic noise  associated  with  them,  are  then  used  to  construct  the  synthetic  sig- 
nals. 

To  study  the  effect  of  the  noise,  the  synthetic  signals  are  con- 
structed as  follows: 
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TABLE  III-  1 

CEPSTRUM  DECOMPOSITION  OF  FOUR  EKAZ  EVENTS 


t(sec): 


m 


o 1 


5 1.0 


FIGURE  111- 1 

WAVEFORMS  OF  SIGNAL,  KA7/081/04N,  AND  CEPSTRUM 
RESOLVED  SIGNALS  ,s  (n)  AND  ,s  p(n-6) 


i 


(in-1) 


xf,(,(n)  - s (n)  - as  (n-n  ) + «f<(n) 

SN  p p ° 

where  s (n)  is  the  resolved  P-phase  from  any  one  of  the  events  listed  in 
P 

Table  111-1,  a and  a are  the  scale  constants,  n^  is  the  delay  time  and 
t(n)  is  the  ncise.  In  equation  (111- 1 ),  -s  (n-nQ)  simulates  the  surface  re- 
flection for  the  ideal  case,  where  the  surface  reflected  pP-phase  is  the  exact 
replica  of  the  P-phase.  Thus,  this  synthetic  mixed  signal  consists  of  two  i- 
dentical  signals  and  some  additive  noise.  Hence,  the  dissimilarity  of  the  P- 
phase  and  the  pP-phase  does  not  enter  into  the  problem  and  we  can  see  the 
sole  effect  of  the  noise.  The  constant  a is  used  to  adjust  the  signal-to-noise 
ratio  (SNR)  such  that  we  can  examine  the  effect  of  various  noise  levels  on  the 
detection  of  the  P-pP  delay  time  and  the  successful  decomposition  of  x^^Jn). 
The  typical  waveforms  and  amplitude  spectra  of  x^fVj(n)  anc*  *(n)  are  S*ven 
in  Figure  111-2  for  SNR  = 00  and  24  dB. 

To  study  the  effect  of  the  dissimilarity,  the  synthetic  signals 
are  formed  as  follows: 

xSN(">  = sl(n)±  s2(n-n0)  (111-2) 

where  s (n)  is  the  resolved  P-phase  and  s (n)  can  be  the  resolved  P-phase 

1 c 

or  pP-phase  from  any  one  of  the  events  listed  in  Table  HI  — 1 • The  positive 
sign  is  used  when  the  resolved  pP-phase  is  assigned  to  s?(n)  and  the  nega- 
tive sign  is  for  the  resolved  P-phase.  Here  i s^n-n^)  simulates  the  im- 
perfect surface  reflection  where  the  surface  reflected  pP-phase  is  not  iden- 
tical to  '.he  direct  P-phase.  The  similarity  coefficient  of  two  time  sequences, 
x(n)  ar.d  y(n)  is  defined  by  (Harley,  1971) 


E(w) 


where  <t>  (0)  is  the  zero  lag  crosscorrelation  value  of  x(n)  and  y(n)  , and 

x y 

,{,  iq)  anc]  (),  (0)  are  the  zero  lag  autocorrelation  values  of  x(n)  and  y(n) 

x-x'  1 yy 

respectively.  The  similarity  coefficient  of  any  two  resolved  phases  is  found 
to  be  in  the  range  of  0.55  - 0.95  (see  subsection  C).  Hen  e,  we  can  show  the 
effect  of  the  dissimilarity  on  the  cepstrun  analyze''  results  by  selecting  ^(n) 
and  s^,(n)  with  varying  similarity  coefficients. 


B.  THE  FFFECT  OF  NOISE 

Four  groups  of  noisy  mixed  signals  were  synthetically  construct- 
ed from  the  resolved  P-phases  of  four  events  listed  in  Table  111-1  (according 
to  the  equation  (111-1)).  In  each  group  there  are  thirty  synthetic  signals  cor- 
responding to  six  noise  levels  and  five  delay  times  (SNR  = ® , 30,  24,  IS,  14, 

b dB  and  n = 4,  5,  6,  7,  8 sampling  units).  Here  the  signal-to-noise  ratio 
o 

is  the  ratio  of  RMS  signal  to  RMS  noise  and  one  sampling  unit  corresponds  to 
0.  1 second.  The  choice  of  the  ranges  for  SNR  and  nQ  is  made  to  reflect  the 
actual  situation  often  encountered  in  the  real  EKAZ  events,  namely: 

• The  NORSAR  recorded  short-period  P-waves  of  EKAZ  events 
usually  have  SNR  between  10  dB  and  35  dB  for  mfa  of  4.  9 to  5.  8. 

• According  to  the  cube  root  scaling  law,  the  P-pP  delay  time  for 

the  contained  explosion  can  range  from  0.3  second  to  0.8  second 

for  m of  4,  5 to  6.  0 (Lane  and  Sun,  1974b). 
b 

• The  shortest  delay  time  that  the  complex  cepstrum  technique 
can  resolve  is  about  0.4  second  (Lane  and  Sun,  1974a,  b). 

The  typical  results  are  presented  in  Figure  HI-3  through  Figure  111-9  for  the 
group  of  syn‘hetic  signals  made  from  the  resolved  P-phase  of  the  event 
KA7./081 /04N.  That  is,  the  synthetic  signals  are 


(n)  = s (n)  - a s (n-fc)  + a ( (n)  i = 1 , 2,  . . . , 6 (111- 3) 

li  1 P 1 P 1 
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CEPSTRA  OF  SYNTHETIC  SIGNAL,  x <n)  WITH  VARYING  SNR 


24  dB 


CEPSTRA  OF  SYNTHETIC  SIGNAL  x (n)  WITH  V, 

DELAY  TIME,  SNR  = 18  dB 


P Similarity  Coefficient 
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30  dB 
24  dB 


18  dB 


12  dB 


5 6 7 8 

n P-pP  Delay  Time 

FIGURE  III- 9 

RITY  COEFFICIENTS  OF  s (n)  AND  s P(n-n  .) 
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where  «.  can  assume  six  values  corresponding  to  six  levels 
i 

of  SNR  and  a is  arbitrarily  assigned  to  be  0.  9 here. 

Figure  111-3  shows  the  cepstra  of  the  signals  with  delay  time 
n = 6 for  six  different  values  of  SNR.  This  figure  illustrates  the  effect  of 
the  noise  on  the  detection  of  the  P-pP  delay  time  for  the  fixed  delay  time. 

In  this  figure,  it  is  noticed  that  the  cepstrak  peaks  due  to  the  P-pP  delay  time 
do  appear  at  n = 6 • k.  k = 1 , 2,  3 as  predicted  in  equation  (11-5).  How  well 
these  cepstral  peaks  can  be  distinguished  from  the  surrounding  cepstra  ob- 
viously depends  on  the  signal-to-noise  ratio.  When  there  is  no  noise  (i.e.  , 

SNR  = oo  ),  the  cepstral  peaks  appear  clearly  at  n = 6,  12,  and  18.  For  SNR 
30  dB,  the  appearance  of  these  cepstral  peaks  remains  more  or  less  the  same 
as  that  for  SNR  = « , except  that  the  surrounding  cepstra  start  growing  and 
interferring  with  the  cepstral  peak  at  n = 18.  The  magnitudes  of  the  surround- 
ing cepstra  keep  increasing  and  the  interference  becomes  more  serious  as  the 
signal-to-noise  rat  o decreases.  Nevertheless,  the  first  two  cepstral  peaks 
can  still  be  distinguished  quite  well  for  SNR  = 24  and  18  dB.  For  SNR  - 12 
and  6 dB,  these  cepstral  peaks  are  buried  in  the  surrounding  cepstra  in  the 
sense  that  they  can  not  be  easily  recogni/.ed  as  in  the  higher  SNR  cases.  Thus, 
strictly  speaking,  for  SNR  = 12  dB,  it  should  be  concluded  that  the  cepstral 
peaks  due  to  the  P-pP  delay  time  cannot  be  detected.  However,  in  most  cases 
analyzed  here,  we  can  detect  these  cepstral  peaks  by  carefully  examining  the 
relative  positions  of  all  suspicious  cepstral  peaks  and  looking  for  the  period- 
icity which  should  exist  for  them  due  to  the  P-pP  delay  time  as  described  by 
the  second  term  in  equation  (11-5). 

Applying  the  shortpass  filter  to  the  cepstra  shown  in  Figure 
111-3  at  n = 6,  the  cepstrum  analysis  decomposes  xSNi(n)  into  the  following 
form 

’‘sNi(n)  = yn,+y>(n-noi>  ■ (UI'4) 
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where  £ (n)  and  £ P(n  - n .)  are  the  cepstrum  resolved  signals  which  are 
p p 01  ° 

expected  to  be  the  estimates  of  s (n)  and  -s  (n-6)  , respectively.  Figure 

p p 

111-4  shows  the  waveforms  of  x (n)  and  'he  cepstrum  resolved  signals.  In 

^5  IN  1 

this  figure,  .s  (n)  is  in  the  solid  line  and  s P(n  - n .)  is  in  the  dashed  line. 

P *~P  ~oi 

It  can  be  seen  that  the  waveform  of  x (n)  remains  very  much  the  same  for 

O IN  1 

SNR  down  to  18  dB.  For  SNR  = « , 30,  24,  and  18  dB,  the  separation  of 

s.  (n)  and  s P(n  - n .)  is  clean  and  the  estimate  of  the  delay  time  n is  6. 
P P oi  x Sii 

which  is  the  actual  delay  time  between  s and  -s  (n-6).  For  SNR  = 12  dB, 

P P 

although  the  separation  is  not  clean,  the  estimate  of  the  delay  time  is  still  6. 
For  SNR  = 6 dB,  it  is  clear  that  the  cepstrum  analysis  fails  to  decompose 
XSNi^  ProPprly-  By  Comparing  equations  ( III- 3 ) and  (III-4),  it  can  be  seen 
that  the  noise  in  xQNI.(n)  is  somehow  distributed  between  £ (n)  and 

pINl  p 

s P(n  - n .).  However,  the  good  resemblance  of  s (n)  and  s (n)  (solid  line 
P oi  "~p  p 

in  Figure  III-2)  and  that  of  s P(n  - n .)  and  -s  (n-6)  (dashed  line  in  Figure 

~~P  ~oi  p 6 

HI- 2 ) are  clearly  observed.  This  indeed  is  quantitatively  reflected  in  the  ap- 
propriate similarity  coefficients  given  in  Table  III- 2 . All  similarity  coeffic- 
ients are  well  above  0.9,  except  those  for  SNR  = 12  dB.  For  SNR  = 1?  dB, 
they  are  still  well  above  0.  8 which  in  general  is  considered  quite  similar.  It 

is  noticed  that  the  similarity  between  £ (n)  and  s (n)  is  better  than  that  be- 

P P 

tween  s P(n  - jn  .)  and  -s  (n-6).  This  might  imply  that  the  noise  is  not 

P Ol  p or/ 

equally  divided  between  the  resolved  signals.  The  second  resolved  signal 

£ P is  noiser. 

P 

Figure  III  - 5 shows  the  complex  cepstra  of  the  signals  with 
SNR  = 18  dB  for  five  different  delay  times.  That  is,  the  synthetic  signals 
are 


*SNi(,') 


s (n)  - s (n  - n .)  + «((n) 
P P oi 


( III  - 5 ) 


where  n = 4,  5,  6,  7,  and  8.  This  figure  shows  the  effect  of  varying  delay 
time  on  the  detection  of  the  P-pP  delay  time  for  the  fixed  noise  level.  There 
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is  a weak  indication  that  the  detection  of  the  eepstral  peaks  due  to  the  P-pP 
delay  time  is  better  for  large  n^.  Nevertheless,  these  eepstral  peaks  can 
be  identified  without  difficulty  for  r^.  down  to  5 and  they  all  appear  at  the 
correct  cepstrum  times,  i.  e.  , n = k%..  k = l,2,...  as  indicated  by  the 
equation  (11-5).  For  n^  4 , the  eepstral  peak  at  n-4  is  concealed  by  the 
cepstrum  of  the  signal  Sp(n).  Hov  ever,  those  at  n = 8 and  12  can  still  be 
identified  with  careful  examination  as  previously  mentioned.  F igure  I1I-6 
gives  the  waveforms  of  ^.(n)  and  the  cepstrum  resolved  signals.  Again, 
the  resolved  signals  closely  resemble  the  original  signals.  The  similarity 
coefficients,  given  in  Table  111-3,  are  all  well  above  0.  9.  There  is  no  de- 
finite relation  between  the  similarity  coefficient  and  the  delay  time.  In  other 
words,  the  longer  delay  time  does  not  necessarily  imply  better  or  worse 
similarity,  one  way  or  the  other.  However,  -he  similarity  between  Sp(n) 

and  s (n)  is  better  than  that  between  s P(n  - n q.)  and  -Sp(n  - n^) 

P ^ 

found  for  the  case  of  fixed  delay  time  and  varying  SNR. 

In  Figures  111-7  through  ill-9,  the  similarity  coefficient  of 
s („)  and  s (n)  , that  of  4pP(n  - nj  and  -Sp(n  - nj.  and  that  of  Sp<"> 
and  s P(n  - n .)  are  plotted  as  functions  of  the  P-pP  delay  time  Op.  using 
signad-to  * noise' ratio  as  the  parameter.  These  figures  display  clearly  the 
relation  amont  the  similarity  coefficient,  the  P-pP  delay  time,  and  the  signal. 

to-noise  ratio. 


as  follows: 


The  major  results  from  the  study  on  the  effect  of  the  noise  are 

The  eepstral  peaks  due  to  the  P-pP  delay  time  can  be  identified 
for  5 ig  nal- to  - noise  ratio  down  to  12  dB  and  the  P-pP  delay  time 
as  short  as  0 4 seconds. 

The  cepstrum  analysis  can  successfully  recover  the  P-phase 
and  pP-phase  for  the  P-pP  delay  time  down  to  0.4  seconds  with 
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TABLE  111-3 

SIMILARITY  COEFFICIENTS  FOR  VARIOUS  P-pP 
DELAY  TIMES  AND  SNR  OF  1 8 dB 


P-pP  Delay  Time 
(Seconds) 

P(s  , s ) 
P P 

PUpp  , -«p> 

% • -V’ 

0.  8 

0.  984 

0.  975 

0.  962 

0.7 

0.  914 

0.917 

0.  879 

0.  6 

0.  965 

0.  950 

0.  956 

0.5 

0.  982 

0.949 

0.  931 

0.4 

0.  977 

0.  945 

0.867 

p(x,y):  similarity  coefficient  of  signals  x(n)  and  y(n) 
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signal-to-noise  ratio  above  18  dB  and  to  0.6  seconds  with 
SNR  = ,'1  dB. 

• For  the  fixed  P-pP  delay  time,  the  detection  of  the  cepstral 
peaks  is  clearly  better  for  high  signal-to-noise  ratio.  The 
similarity  coefficients  increase  significantly  with  the  increas- 
ing signal-to-noise  ratio. 

o For  the  fixed  noise  level,  there  is  only  a weak  indication  that 

the  detection  of  the  cepstral  peaks  is  better  for  large  P-pP  de- 
lay time.  For  high  signal-to-noise  ratio  (SNR  >18  dB),  the 
similarity  coefficients  increase  slightly  as  the  P-pP  delay  time 
increases.  However,  there  is  no  such  relation  observed  for 
SNR  < 12  dB. 

• The  similarity  coefficient  of  and  s^  is  higher  than  that  of 

s P and  -s  . This  implies  that  the  noise  will  seriously  af- 
~~P  P 

feet  the  recovery  of  the  pP-phase. 

C.  THE  EFFECT  OF  DISSIMILARITY 

From  the  eight  cepstrum  resolved  phases  of  the  EKAZ  events 
listed  in  Table  III-  1 , eight  synthetic  signals  are  made  according  to  equation 
(III-2)  for  the  delay  time  of  0.  6 second  (i.  e.  , nQ=  6).  To  facilitate  the  dis- 
cussion below,  a nomenclature  is  established  for  these  synthetic  signals  in 
Table  III-4.  Also  given  in  this  table  are  the  similarity  coefficients  of  s^n) 
and  s2(n).  It  is  noticed  that  these  similarity  coefficients  vary  from  0.55  to 
1.0.  Thus,  the  analysis  of  these  synthetic  signals  simulates  the  situation 
where  the  mixed  signal  consists  of  two  non-identical  signals  with  varying  de- 
gree of  similarity. 

The  cepstra  of  xgD(n)  are  shown  in  Figure  III- 10.  For  the 
signal  SMD-I  which  consists  of  two  identical  signals,  the  cepstral  peaks  due 
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TABLE  111-4 


SYNTHETIC  SIGNAL  CONSISTING  OF  TWO  NON-IDENTICAL 
SIGNALS  WITH  VARYING  DECREE  OF  SIMILARITY 


Signal  I.  D. 

XSD(n)» 

Sj(n) 

s2(n) 

P(sl  , s2) 

SMD-1 

, s (n) 
1 P 

-lSp(n) 

1.0 

SMD-2 

i%(n) 

.S  P(n) 
3 P 

0.933 

SMD-3 

isP(,,) 

2spp(n) 

0.856 

SMD-4 

iVn) 

-3Sp(™) 

0.  844 

SMD-5 

iV"’ 

- s (n) 
4 p 

0.798 

SMD-6 

lspln) 

s P(  n) 
1 P 

0.712 

SMD-7 

lspln) 

-,«(») 
L p 

0.699 

SMD-8 

lsp(n) 

s P(n) 
4 P 

0.  559 

XSD(n)  = Sjt")  + s2(n-6) 


P(Sj  , s^):  similarity  coefficient  of  Sj(n)  and  s^(n) 
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FIGURE  III- 1 0 


CE^STRA  OF  SYNTHETIC  SIGNALS, 
SIMILARITY  BETWEEN  s. 


x (n),  FOR  VARYING 
(n)  AND  sJn) 


to  the  P-pP  delay  time  appear  precisely  at  n - 6k,  k = 1,2,  3 as  Prt" 

dicted  in  equation  (11-10)  when  s^n)  and  s2(n)  are  identical.  However,  as 
soon  as  that  s^n)  and  s^n)  become  non-identical,  even  for  the  signal  SMD-2 
with  a high  similarity  of  0.  933,  the  nice  periodic  occurrance  of  these  cepstral 
peaks  disappears.  Thus,  in  general,  the  cepstral  peaks  due  to  the  P-pP  delay 
time  do  not  exhibit  the  periodicity  as  predicted  by  the  second  term  in  the  equa- 
tion (11-10).  Nevertheless,  a cepstral  peak,  although  not  as  sharp  as  theoret- 
ically expected,  is  observed  at  n=6  where  the  first  cepstral  peak  should  appear. 
For  a similarity  coefficient  greater  than  0.7,  the  identification  of  this  cepstral 
peak  does  not  become  more  difficult  as  the  similarity  coefficient  decreases. 

Applying  the  shortpass  filter  to  the  cepstra  at  n = 6 , the  cep- 
strum  analysis  decomposes  these  synthetic  signals  into  the  following  form 


XSD(n)  = 'l(n)  + f°r  SMD_i’  { = 1,8 


(111-6) 


The  waveforms  of  these  signals  are  given  in  Figure  111-11.  It  can  be  seen  that 
the  resemblance  between  s^n)  and  s^n)  . and  that  between  ^(n)  and  s.,(n) 
are  good.  This  is  clearly  shown  by  their  similarity  coefficients  given  in  Table 
111-5.  By  comparing  the  positions  of  their  corresponding  amplitude  peaks,  the 
delay  time  between  s^n)  and  s^n)  is  estimated  to  be  0.6  seconds  as  indi- 
cated in  equation  (111-6).  This  is  exactly  equal  to  that  between  s^n)  and  s^n). 
Therefore,  cepstrum  analysis  successfully  decomposes  xgE>(n)  into  its  con- 
stituent parts. 

The  major  results  from  the  study  on  the  effect  of  dissimilarity 
can  be  summarized  as  follows: 

• The  dissimilarity  between  s^n)  and  s^n)  will  destroy  the 

periodicity  of  the  cepstral  peaks  due  to  the  P-pP  delay  time. 
However,  the  first  cepstral  peak  does  occur  at  the  right  cep- 
strum time  which  is  equal  to  the  P-pP  delay  time. 
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TABLE  III  - 5 


SIMILARITY  COEFFICIENTS  FOR  VARYING  SIMILARITY 
BETWEEN  s j (n)  AND  s2(n) 


Signal  I.  D.  p^  , s2)  P^  , ^ P(s{  , s ^ p(s,2  , s2) 


SMD-1 


SMD-2 


SMD-3 


SMD-4 


SMD-5 


SMD-6 


SMD-7 


SMD-8 


0.  933 


0.  856 


0.  844 


0.798 


0.712 


0.699 


0.  559 


0.  938 


0.884 


0.  870 


0.781 


0.  723 


0.  735 
0.  564 


0.991 


0.997 


0.992 


0.995 


0.918 


0.891 


0.871 


0.  991 


0.  995 


0.  992 


0.  998 


0.890 


0.  826 


0.833 


p(x,  y):  similarity  coefficient  of  x(n)  and  y(n) 


• The  identification  of  the  first  cepstral  peak  is  independent  of 
the  similarity  between  Sj(n)  and  s^n).  There  is  only  a very 
weak  indication  that  lower  similarity  may  result  in  poor  iden- 
tification. 

• The  cepstrum  resolved  signals,  .s^n)  and  s^n)  , agree  very 
well  with  the  original  signals,  s^(n)  and  s^(n).  Their  sim- 
ilarity coefficients  are  well  above  0.  8. 

• The  identification  of  the  first  cepstral  peak  and  the  successful 
decomposition  of  Xg^(n)  by  apply*ns  the  shortpass  filter  at  the 
occurrance  of  this  cepstral  peak  will  constitute  the  detection  of 
the  P-pP  delay  time.  In  this  sense,  cepstrum  analysis  can  de- 
tect the  P-pP  delay  time  and  successfully  decompose  xSD(n) 
for  a similarity  coefficient  of  s^(n)  and  s^n)  aS  l°w  as  0.55. 

By  comparison  of  the  results  obtained  here  and  those  in  the  Subsection  1II-E, 
it  is  reasonable  to  say  that  the  dissimilarity  will  affect  the  detection  of  the 
cepstral  peaks  due  to  the  P-pP  delav  time  more  seriously  than  does  the  noise. 
In  general,  for  a mixed  signal  consisting  of  two  non-identical  signals,  we 
should  not  expect  to  find  .ne  periodicity  of  the  cepstral  peaks  due  to  the  P-pP 
delay  time.  Father,  we  should  look  for  the  suspicious  cepstral  peaks  and 
apply  the  shortpass  filters  accordingly.  The  correct  estimate  of  the  P-pP  de- 
lay time  can  only  be  obtained  through  the  successful  cepstrum  decomposition 
of  the  mixed  signal  Xg^(n). 
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d,  the  effect  of  prefiltering 

It  has  been  shown  in  Subsection  B that,  to  some  degree,  the 
noise  affects  the  detection  of  the  cepstrai  peaks  due  to  the  P-pP  delay  time. 

The  purpose  of  this  experiment  is  to  find  out  if  the  reduction  of  the  noise  by 
filtering  before  cepstrum  analysis  will  improve  detection.  The  experiment 
is  carried  out  on  the  same  synthetic  signals  as  those  used  in  Subsection  B. 

From  the  amplitude  spectra  of  the  noise  and  the  synthetic  signal  (SNR  = « ) 
shown  in  Figure  III-2a.  it  is  noticed  that  the  noise  spectrum  is  confined  be- 
low 2.  5 Hz  and  the  signal  has  significant  spectral  content  from  0.  5 Hz  to  the 
folding  frequency  which  is  5 Hz.  Thus,  it  is  clear  that  there  is  some  over- 
lapping of  the  noise  and  the  signal  spectra,  although  their  main  spectra  are 
widely  separated,  0-0.75  Hz  versus  1.25-3  Hz. 

The  filter  which  is  employed  here  is  the  zero  phase  bandpass 
filter  with  cosine  tapering  at  both  ends,  as  shown  schematically  in  Figure 
III- 12.  The  cut-off  and  corner  frequencies  at  the  low  and  high  ends  of  the 
filters  are  represented  as  f^,  f L*  and  f^  fj_j0*  respectively, 
groups  of  filters  are  used  here.  Their  characteristic  frequencies,  f^.  t ^ 
f , and  f are  given  in  Table  III-6.  The  first  group  is  the  highpass  filter 
with  varying  low  frequency  cut-off.  This  is  intended  to  show  the  effect  of  re- 
ducing the  noise  which  mostly  is  in  the  low  frequency.  The  second  group  is 
the  lowpass  filter  with  varying  high  frequency  cut-off.  This  is  to  show  the 
effect  of  the  minor  change  in  the  signal  spectral  shape.  The  last  group  is 
the  bandpass  filter  . show  the  combined  effect  of  the  noise  reduction  and  the 
minor  alteration  of  the  signal  spectral  shape. 

The  typical  results  are  presented  here  for  the  synthetic  signal 
given  by  equation  (HI-3)  with  SNR  = 18  dB.  Figured  111-13  through  111-19  show 
the  results  using  the  first  group  of  filters.  For  eachftgure.  the  cepstrum 

and  the  amplitude  spectrum  of  the  filtered  signal.  *SN(")  • are  glven  £irSt' 
followed  by  the  waveforms  of  the  filtered  signal  and  the  cepstrum  resolved 
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FIGURE  III- 1 2 

RESPONSE  OF  THE  BANDPASS  FILTER 


TABLE  ill -6 

SPECIFICATIONS  OF  THE  BANDPASS  FILTERS 
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FIGURE  III- 13 

CEPSTRUM  ANALYZED  RESULTS:  PREFILTERING  WITH  FILTER 


CEPSTRUM  ANALYZED  RESULTS:  PREFiLTERING  WITH  FILTER  I-Z 


FIGURE  TI-16 

CEPSTRUM  ANALYZED  RESULTS:  ^RE  FILTERING  WITH  FILTER  1-4 


FIGURE  III—  1 7 

CEPSTRUM  ANALYZED  RESULTS:  PREFILTERING  WITH  FILTER  I- 


signals,  sF(n)  and  s*P(n-nJ.  For  comparison  purposes,  results  without 
h ~p  ~P  ° 

prefiltering  are  presented  in  Figure  111-20.  It  is  clear  that  no  visible  change 
in  the  waveform  of  the  filtered  signal  is  observed  as  compared  to  that  of  the 
unfiltered  signal,  except  for  the  last  filter  whose  an<^  into  the 

main  part  of  the  signal  spectrum.  However,  from  their  amplitude  spectra, 
we  can  see  that  the  low  frequency  noise  is  reduced  and  the  shape  of  the  low 
frequency  spectrum  is  deformed.  Nevertheless,  the  general  shape  of  the 
spectrum,  especially  the  occurrances  of  the  spectral  maxima  and  minima, 
remains  very  much  the  same  for  a low  frequency  cut-off  up  to  1.0  Hz.  The 
cepstral  peaks  due  to  the  P-pP  delay  time  appear  at  the  right  cepstrum  times 
(n  = 6k,  k = 1,2,...)  and  can  be  identified.  It  is  hard  to  judge  if  the  detec- 
tion of  these  cepstral  peaks  is  better  for  the  filtered  signals.  It  is  fair  to  say 
that  the  prefiltering  with  the  low  frequency  cut-off  up  to  1.0  Hz  neither  im- 
proves nor  downgrades  the  detection.  Also,  the  cepstrum  resolved  signals 
of  the  filtered  and  the  unfiltered  signals  are  very  much  the  same.  When  the 
lower  end  of  the  filter  is  moved  up  to  1.5  Hz,  the  waveform  and  the  low  fre- 
quency spectrum  of  the  filtered  signal  differ  significantly  from  those  of  the 
unfiltered  signal.  Although  the  first  cepstral  peak  still  appears  at  the  right 
cepstrum  time,  the  identification  of  the  second  cepstral  peak,  which  is  sup- 
posed to  be  at  n = 12,  can  no  longer  be  made. 

Figures  111-21  through  111-25  present  similar  results  using  the 
second  group  of  filters.  It  is  noticed  that  there  is  a long,  low  level  high  fre- 
quency tail  (3.5  Hz  - 5 Hz)  in  the  amplitude  spectra  of  the  unfiltered  signal. 
Again,  the  waveforms  of  the  filtered  signals  are  all  very  similar  to  that  of  the 
unfiltered  signals.  However,  when  the  high  frequency  cut-off  is  below  4.  5 Hz 
some  of  this  high  frequency  tail  is  eliminated,  and  the  second  cepstral  peak 
due  to  the  P-pP  delay  time  is  no  longer  there.  Thus,  even  a very  minor 
change  of  the  signal  spectrum  will  seriously  affect  the  detection  of  the  P-pP 

delay  time. 
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CEPSTRUM  ANALYZED  RESULTS:  NO  PREFILTERING 


CEPSTRUM  ANALYZED  RESULTS:  ^RE FILTERING  WITH  FILTER  II-2 


FIGURE  III— 2 3 

CEPSTRUM  ANALYZED  RESULTS:  °RE  FILTERING  WITH  FILTER  II- 


FIGURE  III — 2 5 

CEPSTRUM  ANALYZED  RESULTS:  PREFILTERING  WITH  FILTER  II- 


The  result  (or  the  last  Broup  is  give-  in  Figure  W-2<>.  Th'5 
U.,  filter  only  passes  the  main  part  of  the  signal  spectrum.  This  simulates 
the  normal  operation  of  the  bandpass  filter  to  improve  the  s.gnal- to- no.se 
ratio.  The  waveform  of  the  filtered  signal  again  is  very  similar  to  that  of  t e 
unfilte red  signal,  but  the  cepstrum  analyzed  results  are  worse  than  those  for 
the  second  group.  The  cepstral  peahs  due  to  the  P-pP  delay  time  cannot  be 
identified  and  the  cepstrum  decomposition  is  incorrect. 

The  major  results  from  this  study  on  the  effect  of  prefiltering 

Prefiltering  with  a bandpass  filter  does  not  affect  the  cepstrum 
analyzed  results  as  long  as  the  passband  of  the  filter  covers 
the  entire  signal  spectrum. 

The  reduction  of  the  noise  by  prefiltering  does  not  compensate 
the  effect  of  minor  distortion  of  the  amplitude  spectrum  (as- 
sociated  with  the  prefiltering)  on  the  detection  of  the  P-pP  de- 
lay time. 


are  as  follows: 
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CEPSTRUM  ANALYZED  RESULTS:  PREFILTERING  WITH  FILTER 
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SECTION  IV 

EXPERIMENTAL  RESULTS:  REAL  SIGNALS 


I 
I 
I 
I 

| In  this  section  cepstrum  analysis  is  applied  to  several  short- 

period  P-waves  from  EKAZ  events  recorded  at  NORSAR.  Table  IV  - 1 pre- 
sents source  information  for  these  events,  including  the  calculated  signal- 
to-noi3e  ratio.  The  signal-to-noise  ratios  for  these  events  range  from  ’0 
dB  to  40  dB. 

Before  applying  the  cepstrum  analysis  to  these  signals,  it  is 
desirable  to  set  up  some  criteria  for  judging  the  results.  To  facilitate  the 

statements  of  these  criteria,  let  us  assume  that  the  signal  x(n)  is  cepstrum 

i 

decomposed  into  two  resolved  signals  as  follows: 


I 

I 


I 

I 

I 

I 

I 

I 

1 


x(n)  = .s^n)  + s2(n  - nQ)  . 

A real  signal  will  be  said  to  be  successfully  decomposed  by  the  cepstrum  anal- 
ysis if: 

m s (n)  and  s (n  - n ) resemble  a relistic  teleseismic  short- 

' ' — i ' ' ~-2  o 

period  P-phase  of  a single  underground  explosion  (i.  e.  , if  they 
are  simple  and  have  reasonably  finite  trace  lengths).  Also,  the 
first  few  significant  amplitude  peaks  of  s^n  “ s^ould 

delayed  by  n samples  with  respect  to  the  corresponding  peaks 
of  js  (n);  and  if  _s  is  the  surface  reflection  of  £ , they  should 

be  of  opposite  sign. 

(2)  s (n)«0  for  n > Nj  , where  N^  is  a reasonably  finite  integer 

(for  the  EKAZ  events  recorded  at  NORSAR,  Nj  = 1.  5-2.0  sec- 
onds). 
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TABLE  IV - 1 

PRESUMED  UNDERGROUND  EXPLOSIONS  in  EASTERN  KAZAKH 


Event 
I.  D. 

Location 

mb 

SNR 

Date 

Latitude 

Longitude 

EKZ/  159/01N 

06/07/72 

49.  8 

78.2 

5.5 

39 

KA7/081/04N 

03/22/71 

49.7 

78.2 

5.8 

34 

EK7  /239/03N 

08/26/72 

SO.  0 

77.  8 

5.  5 

33 

KA7./170/04N 

06/19/71 

50.  0 

79.  1 

5.  5 

30 

EK7  /333/06N 

11/29/71 

49.  8 

78.  1 

5.  5 

28 

KAZ/282/06N 

10/09/71 

50.  0 

77.7 

5.4 

22 

EK7  /246/03N 

09/02/72 

50.  0 

77.  7 

5.  1 

14 

EK7/363/04N 

12/28/72 

50.  0 

78.  0 

4.5 

. 

10 

i 


I 

I 

I 

I 

I 
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I 

I 

I 
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(3)  s2(n  - n ) =0  for  n < nQ  . 

Figures  1V-1  through  1V-8  present  the  results  of  cepstrum 
analyzing  these  signals.  For  each  figure,  the  waveform  and  the  cepstrum  of 
the  signal  are  displayed,  together  with  the  cepstrum  resolved  P-phase  (solid 
line)  and  the  pP-phase  (dashed  line)  By  examining  these  cepstra,  even  for 
those  high  SNR  events,  we  cannot  find  the  periodic  occurrances  of  the  cep- 
stral  peaks  due  to  the  P-pP  delay  time  as  that  discussed  in  Subsection  III-B. 
This  is  thought  to  be  reasonable,  since  in  general  the  surface  reflected  pP- 
phase  is  not  the  exact  replica  of  the  direct  P-phase.  However,  it  is  noticed 
that,  except  for  two  low  SNR  events,  the  appearances  of  these  cepstra  are 
similar  to  those  demonstrated  in  Subsection  1II-C,  where  the  synthetic  signal 
consisted  of  two  non-identical  signals.  Thus,  for  the  real  signals,  to  obtain 
the  estimate  of  the  possible  P-pP  delay  time,  the  successful  decomposition 
of  these  signals  must  be  achieved  first  as  mentioned  in  Subsection  III-C. 

Since  the  cepstrum  does  not  show  the  definite  peaks  due  to  the 
possible  P-pP  delay  time,  the  most  straightforward  way  to  filter  the  cepstrum 
is  to  apply  the  shortpass  filter  at  every  cepstral  peak  which  is  suspected  to  be 
the  cepstral  peak  due  to  the  P-pP  delay  time.  Then  the  corresponding  op- 
strum  resolved  signals  are  examined  to  see  if  the  signal  is  successfully  de- 
composed. As  an  example,  consider  the  event  EKZ/333/06N  as  shown  in 
Figure  IV -5.  The  cepstral  peaks  at  n = 5,  6,  7,  8,  and  10  are  decided  to  be 
the  suspicious  peaks.  The  cepstrum  decomposition  by  applying  the  shortpass 
filter  at  n = 7 is  shown  in  Figure  IV-5a  and  is  judged  unsuccessful.  The 
same  conclusion  is  made  for  n = 8 and  10.  The  cepstrum  resolved  signals 
shown  in  Figure  IV-5  are  obtained  by  applying  the  shortpass  filter  at  n = 6. 
The  similar  successful  decomposition  is  also  achieved  by  filtering  at  n = 5. 
The  P-pP  delay  time  is  estimated  to  be  0.6  seconds. 
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CEPSTRUM  ANALYZED  RESULTS 


FIGURE  IV -6 

CEPSTRUM  ANALYZED  RESULTS:  KAZ/282/06N 


The  final  cepstrum  decompositions  of  these  signals  are  expres- 
sed in  the  form  like  equation  (1II-6)  and  are  given  in  Table  1V-2.  Also  shown 
are  the  similarity  coefficient  and  the  amplitude  ratio  of  the  resolved  signals. 

It  is  noticed  that  these  similarity  coefficients  range  from  0.7  to  0.9  which 
are  compatible  with  those  of  the  synthetic  signals  used  in  Subsection  III-C  . 
However,  the  trend  relation  of  increasing  similarity  with  increasing  signal- 
to-noise  ratio,  which  has  been  found  in  Subsection  III- B for  the  synthetic  sig- 
nal with  identical  P-phase  and  pP-phase,  is  not  observed  here. 
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TABLE  IV-2 

PINAL  CEDSTRUM  DECOMPOSITIONS  OF  EKAZ  EVENTS 


Eve  nt 
I.  D. 

Cepstrum  Decomposition 

(1) 

a 

(2) 

P 

EKZ/159/01N 

x(n)  = s (n)  + s P(n-7) 
P P 

0.  87 

0.  864 

KAZ/081/04N 

x(n)  = s (n)  + s P{n-6) 
P P 

0.  74 

0.  712 

EKZ/239/03N 

x(n)  = s (n)  + s P(n-5) 
P P 

0.  60 

0.  690 

KAZ/170/04N 

x(n)  = s (n)  + s P(n-6) 
P P 

0.  83 

0.796 

EKZ/  333/06N 

x(n)  = s (n)  + s P(n-6) 
P P 

I.  05 

0.  906 

KAZ/282/06N 

x(n)  = s (n)  + s P(n-5) 
P P 

0.  85 

0.  874 

EKZ/246/08N 

x(n)  = s (n)  + s P(n-6) 
P P 

0.  65 

0.749 

EKZ/363/04N 

x(n)  = s (n)  + s P(n-5) 
P P 

1.  03 

0.  730 

(1)  a : peak-to-peak  amplitude  ratio  of  s P(n)  and  s (n) 

P P 


(2)  p : similarity  coefficient  of  s (n)  and  s P(n) 

P P 


f 
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SECTION  V 
CONC  LUSIONS 


The  noise,  the  dissimilarity  between  the  P-phase  and  the  pP- 
phase,  and  the  prefiltering  are  found  to  produce  various  unfavorable  effects 
on  the  detection  of  the  P-pP  delay  time  and  the  successful  cepstrum  analysis 
of  the  mixed  P-  and  pP-wave.  The  dissimilarity  will  affect  the  detection 
more  seriously  than  the  noise  does.  The  increase  of  the  noise  level  will  only 
gradually  obscure  the  cepstral  peaks  due  to  the  P-pP  delay  time,  while  the 
dissimilarity  completely  destroys  the  periodic  occurrances  of  these  cepstral 
peaks.  For  the  identical  P-phase  and  the  pP-phase,  the  cepstrum  analysis 
can  detect  the  P-pP  delay  time  as  short  as  0.4  second  and  successfully  re- 
cover them  for  the  signal-to-noise  ratio  as  low  as  12  dB.  For  the  non- 
identical P-phase  and  pP-phase,  the  cepstrum  analysis  can  still  achieve  the 
successful  decomposition  for  the  similarity  coefficient  about  0.  55.  This,  in 
turn,  can  be  used  to  estimate  the  P-pP  delay  time. 

In  reality,  the  pP-phase  and  the  P-phase  are  not  identical,  al- 
though they  can  be  expected  to  be  quite  similar.  Therefore,  for  the  real  mix- 
ed P-  and  pP-wave,  we  should  not  expect  to  find  the  periodicity  of  the  cepstral 
peaks  due  to  the  P-pP  delay  time.  Rather,  we  should  look  for  the  suspicious 
cepstral  peaks  and  filter  the  cepstrum  accordingly  for  the  possible  decom- 
position. The  correct  estimate  of  the  P-pP  delay  time  can  only  be  obtained 
through  the  successful  cepstrum  decomposition  of  the  mixed  signal. 

Prefilteriug,  in  general,  does  not  affect  the  cepstrum  analyzed 
results  as  long  as  the  passband  of  the  filter  covers  the  entire  signal  spectrum. 
However,  it  is  suggested  that  the  prefiltering  should  be  avoided  if  the  signal- 
to-noise  ratio  is  not  unreasonably  low. 
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The  effects  of  the  noise,  the  dissimilarity,  and  the  prefiltering  on  the 
detection  of  the  P-pP  delay  time  and  the  successful  decomposition  of  the  mix- 
ed P-  and  -pP  wave  by  employing  the  cepstrum  analysis  arc  studied  here. 

These  effects  have  been  investigated  by  extracting  the  direct  P-phascs  and  the 
surface  reflected  pP-phas us  from  several  presumed  underground  explosions 
in  eastern  Kazakh,  and  remixing  them  at  various  time  delays  with  the  addition 
of  various  levels  of  the  realistic  seismic  noise.  The  quality  of  the  cepstrum 
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analyzed  'osuUs  as  a function  of  the  signal-, o-noisc  ratio,  the  delay  time, 
and  l je  dissimilarity  between  the  P.phase  and  Ine  pP-phase  have  been  ob- 
tamed  for  a number  of  such  events.  In  general,  it  is  found  that  tne  c,  r- 
strum  analysts  can  detect  the  P-pP  delay  time  as  short  as  0.1  second  and 
successfully  recover  the  P-  and  the  pP-phases  for  the  signal- 'o-noise 
rat.o  down  to  about  12  dB.  When  the  similarity  coefficient  between  the 
. ” and  the  Pp-Phases  is  less  than  1.0,  the  detection  of  the  P-pP  delay 
time  becomes  not  so  obvious  as  it  is  when  they  are  identical.  The  dis- 
simUar^y  between  the  P-phase  and  the  pP-phase  will  affect  the  detection 
of  the  P-pP  delay  time  more  seriously  than  the  noise  in  the  cepstrum 
analysis.  Prefiltering  of  the  noisy  mixed  P-  and  -pP  wave  with  the 
standard  bandpass  filter  will  not  affect  the  cepstrum  analyzed  results  as 
ong  as  the  filter  band  covers  the  entire  signal  spectrum;  although  it  seems 
that  no  prefiltenng  will  offer  better  detection  of  the  P-pP  delay  time. 
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